Computer Assembly of Cluster-Forming Amphiphilic Dendrimers 
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Recent theoretical studies have predicted a new clustering mechanism for soft matter particles 
that interact via a certain kind of purely repulsive, bounded potentials. At sufficiently high densities, 
clusters of overlapping particles are formed in the fluid, which upon further compression crystallize 
into cubic lattices with density-independent lattice constants. In this work we show that amphiphilic 
dendrimers are suitable colloids for the experimental realization of this phenomenon. Thereby, we 
pave the way for the synthesis of such macromolecules, which form the basis for a novel class of 
materials with unusual properties. 

PACS numbers: 82.70.Dd, 61.20.Ja, 64.75.Yz, 81.05. Zx 
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Self-assembly is the key-expression that circumscribes 
the incredibly rich wealth of ordered phases encoun- 
tered in soft matter systems. Starting from the face- 
center cubic arrangement assumed by colloidal spheres 
at high concentrations P, [2j 01 E^nd low-symmetry crys- 
tals formed by soft spheres [4|, they extend over to a 
variety of alloys seen for charged colloidal mixtures @ 
and to the gyroid phases assembled in block copolymer 
solutions [l,13i 9 • Recent theoretical and computational 
advances have predicted a novel form of self-assembly in 
soft matter, i.e., the formation of stable clusters 0, 
encountered for purely repulsive, bounded effective po- 
tentials <i>off(r). These clusters crystallize into cubic lat- 
tices at sufficiently high densities and all temperatures 
0i [lot ■ This phenomenon bears significant consequences 
both from the fundamental point of view [l^ and from 
the aspect of the properties of the ensuing materials, e.g., 
their diffusion and relaxation dynamics [Tl| . Though 
thoroughly understood at the level of effective potentials, 
the phenomenon begs the question: What kinds of parti- 
cles display the class of effective interactions giving rise 
to this phase behavior? In this contribution we demon- 
strate that relatively simple macromolecules can be de- 
signed to achieve this aim, taking advantage of the great 
fiexibility offered by soft matter systems to manufacture 
new materials. 

The underlying theoretical concepts can be clearly 
stated as follows. Repulsion-induced aggregation which 
leads to ordered cluster phases requires that the Fourier 
transform (FT) "I>cff(fc) of ^csir) has negative parts for 
some values of the wave- vector fc. If, however j_$pg(fc) > 
for all fc, reentrant melting occurs instead [l2|- A suf- 
ficient condition for the former is that ^'^g{r = 0) > 
[lot . Useful archetypes of bounded interactions ^, To'l are 
the generalized exponential models of index n (GEM-n), 
where <l?„(r) = eexp[— (r/cr)"], with e and a being some 
energy- and length-scales: here, for n > 2 clustering takes 
place, whereas for n <2 reentrant melting occurs [isj . 

Searching for realizations of clustering-type potentials, 
we concentrate on dendrimers, a choice motivated by 




FIG. 1: (Color online) Schematic representation of the ef- 
fective potentials of amphiphilic dendrimers: interactions 
between the different regions (solvophobic core - yellow, 
solvophilic shell - blue) of amphiphilic dendrimers lead to in- 
teractions that are steeper (red line) than the threshold of 
clustering (Gaussian, black line). 



their outstanding properties. They are characterized 
by a high degree of monodispersity and a well-defined, 
highly branched internal structure; efficient dendrimer 
assembly has been boosted by recent progress in syn- 
thetic techniques [l3] ■ Fundamentally, they serve as tun- 
able soft colloids that allow for control of their effective 
interactions via changes in chemicjil composition, bond 
length and generation number [isl . Il6| . For athermal 
dendrimers, $eff('') has a Gaussian shape (n = 2 in 
the GEM-n family), hinting thus to the reentrant melt- 
ing scenario. Closely linked with these findings is the 
dense-core structure of these dendrimers, arising from 
back- folding of the terminal groups [la, [i3] ■ 

Following the ideas of material design, we modify the 
architecture of athermal, fiexible dendrimers along well- 
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defined strategies. Since the Gaussian effective interac- 
tion is at the threshold to clustering, we require the fol- 
lowing changes in ^cs{r) to achieve clustering behavior: 
a flatter core region, such as those of the GEM-n po- 
tentials with n > 2, which - compared to a Gaussian - 
also display a steeper decay of the repulsion for larger 
separations. Alternatively, a positive effective interac- 
tion with a local minimum at r = also leads to oscil- 
lations in $eff(g), as we have $"fr('' = 0) > in that 
case. To realize this goal, we aim for a more open struc- 
ture and stronger segregation between outer and inner 
particles by assembling amphiphilic dcndrimcrs built up 
from a solvophilic shell and solvophobic core particles. 
In Fig. [ll we qualitatively illustrate how our modifica- 
tions lead in the correct direction: as the macromolecules 
start to overlap, the solvophilic and thus mutually repul- 
sive shells cause a steeply increasing potential wall. This 
effect is reinforced upon further decreasing the distance 
since core and shell repel each other due to their different 
nature. Eventually, the attractive core regions overlap 
and slow down further growth of the repulsion, leading 
to a rather flat region or even a local minimum in ^csif) 
at small distances. 

In an effort to realize these ideas we developed a com- 
puter model of second generation amphiphilic den- 
drimers [l^ where the end-groups form the solvophilic 
shell (index S) and all other monomers the solvophobic 
core (index C). The bonds between monomers are mod- 
eled by the finitely extensible nonlinear elastic (FENE) 
potential 



10 



/3*™M = -^..^' log 




, (1) 



with fiv = CC, CS, which restricts the bond-length to be 
. K^^ is the spring-constant and 



in between Z^'' and 

Ru. = e.f" - /°.„ with /" = {l]^^^^ + /;°j")/2 being the 



equilibrium bond length. All other interactions between 
two monomers separated by distance r are modeled by 
the Morse potential 



- fJ.V 



1 
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(2) 



with ^1/ — CC, CS, SS, which is characterized by a re- 
pulsive core at short and an attractive tail at long dis- 
tances whose depth and range are parametrized by s^i, 
and a^y, respectively. The d^jy are the monomer diam- 
eters. All potential parameters of the dendrimers dis- 
cussed in this work are summarized in Table [H 

We calculated the monomer density profiles for the 
core and the shell particles, pci^) and ps{i'), in stan- 
dard Monte Carlo (MC) simulations of an isolated den- 
drimer. Representative results pertaining to two (I?i, 
D2) out of seven model dendrimers simulated are shown 
in Fig. [21 They demonstrate that both the core- and 
the shell-particle distributions are of Gaussian shape, 
with the fitting parameters given in the respective pan- 
els. In striking contrast to athermal dendrimers (see Fig. 
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FIG. 2: (Color online) Monomer density profiles of the core 
(red) and the shell (green) region as obtained by MC simu- 
lations (symbols) and from fits to Gaussians (lines) for the 
amphiphilic dendrimers Di (left panel) and D2 (right panel). 
The interaction parameters are given in Table I. In the in- 
set, we show the same for an athermal dendrimer of the same 
generation. 
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CS 


0.014 


19.2 


1.50 (Di) 
1.25 {D2) 


SS 


0.014 


19.2 


2.0 (Di) 
1.5 (D2) 



TABLE I: Potential parameters of the dendrimers considered 
in this study [of. Eqs. JT} and labeled Di and D2. ZZ 
refers to the two central monomers. 



[21 inset) , these density profiles are - due to amphiphilic- 
ity - spatially segregated: the distribution of the core 
monomers has its maximum close to the origin while the 
profile of the shell particles is centered between 0.85 i?g 
and Rg (for all seven dendrimers simulated), with Rg 
being the dendrimer's radius of gyration. 

Next, we performed MC simulations between two in- 
teracting dendrimers, averaging over the degrees of free- 
dom of the constituent monomers to determine the ef- 
fective potential, $off(?'). Allowing two dendrimers to 
interact freely, their effective potential can be deter- 
mined from the centers-of-mass (c.o.m.) correlation tunc- 
tion G(r), given by G(|Ri - R2I) = (ei(Ri)^2(R2)). 
Here, gi(Ri) — S^Ri — Si) is the density operator for 
the c.o.m. of dendrimer i, and the average (■ ■ • ) is taken 
over all c.o.m. position vectors Si {i = 1,2). ^cs{r) is 
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FIG. 3: (Color online) Simulation snapshot of two interacting 
amphiphilic dendrimers, both showing a dense shell confor- 
mation. 



then given by l3<^cfiir) = - \n[G{r)], with (3 = l/fceT. 
Since the repulsion between the dendrimers is expected 
to be strong at short distances, this scheme will pro- 
vide only poor statistics for small separations. To cope 
with this problem, we use non-Boltzmann sampling [20| 
where we divide the r-range into m = 15 windows of 
width Ar," = y*""". . ( j + 1) and assume for each window 
J = 0, . . . , 



m(m+l) 

— 1 an umbrella potential Wj (r) 





oo else 



Sr < r < ri 



Ar 



(3) 



Here, rp = and else rj = J2i=o Further, rmax = 

5Rg and dr is chosen to guarantee for slightly overlapping 
windows (at the edges, Sr = 0.) For each of these win- 
dows, simulations of 2 x 10* MC sweeps are performed 
and Gj (r) is determined to within an additive constant 
due to normalization. We obtain $offj ('') in each win- 
dow j as —fee Tin Gj{r) and these results are merged to 
form a continuous curve which is finally normalized by 

setting /3$eff(r- = Tmax) = 0. 

A typical simulation snapshot of two interacting den- 
drimers is shown in Fig. [31 Results for the effective 
interactions of the two dendrimers are summarized in 
Fig. m $cff('^) indeed shows a steep increase as the 
macromolecules approach and eventually becomes rather 
flat [Fig. 4(a)], or even exhibits a locally attractive dip 
[Fig. 4(c)] for smaller distances. We fitted ^cii{r) of den- 
drimer Di by a GEM-n potential, finding that <l>„(r) 
with an index n = 3.1 approximates it with high ac- 
curacy. The effective potential of dendrimer D2 can be 
fitted to a double-Gauss interaction of the form $off ('') — 
ei exp[-(r/cri)^] -62 exp[-(r/(T2)^], where ei = 23.6/^37", 
£2 = 22.5fcBT, CTi = 1.117i?g and 0-2 = 1.059i?g. Both ef- 
fective potentials lead to clustering, since $"ff (r = 0) > 0. 

To provide a semi-quantitative theoretical background 




FIG. 4: (Color online) The effective potentials ^efi(r) of typ- 
ical amphiphilic dendrimers [(a), (c)] and their FTs ^^si^q) 
[(b), (d)], showing negative parts. The blue dashed line de- 
notes the theoretical result (see text), and the green lines are 
fits to the simulation data, (a) and (b) pertain to dendrimer 
Di, (c) and (d) to D2- The insets in (a) and (c) feature 
simulation snapshots of the respective dendrimers. 



to these simulation results, we reconsider our amphiphilic 
dendrimers within a suitably modified Flory theory. The 
original idea of this concept [21|] is based on the simpli- 
fying assumptions that the spherosymmetric monomer 
densities of isolated athermal dendrimers around the 
c.o.m., p{r), do not change upon close interaction. Here, 
p{r) = ~ ^j)) and Tj denotes the position vector 

of monomer j with respect to the center of mass. Thus, 
considering two such dendrimers at separation R and as- 
suming that their profiles are not distorted by their mu- 
tual presence, the effective interaction between the two 
dendrimers takes the form 

$eff(|R|)= |yp(|ri|)p(|r2-R|)«(|ri-r2|)dridr2, (4) 

where w(|ri — r2|) is the monomer-monomer interaction. 
For the latter, a contact interaction weighted by the sec- 
ond virial coefficient vq of the monomer-monomer inter- 
action is introduced, Pv{\ri — r2|) = wo^(|ri — r2|)- For 
athermal dendrimers, vq > 0. The FT of ^cs{t) given in 
Eq. (dl then reads in this case /3<loff (fc) = VQp'^{k) > Vfc. 

Generalizing this model to amphiphilic dendrimers, 
we treat the core and the shell profiles separately, in- 
troducing three different excluded volume parameters, 
i^cCi i^cs^ a-nd wss given by the second virial coefficients 
of the underlying monomer-monomer interactions. Pro- 
ceeding along similar lines as above leads to the FT of 
the effective interaction between amphiphilic dendrimers 
as P^es{k) = Y.ti,v'"i^i'hO^)pAk) where pp(fc) is the 
FT of p^{r) and i^,iy = C, S. Core solvophobicity im- 
plies vqc < 0, whereas wcSjI'SS > 0. Consequently, 
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$cff(fc) can also display negative components. The val- 
ues of the second virial coefficients for our dendrimers are 
vcc = -1.90, vss = 28.40 and vcs = 11-31 (fi), and 
vcc — —1.90, vss — 11.31 and vcs — 6.25 {D2), mea- 
sured in d^Q. Based on the simulation results (cf. Fig. [2]), 
we model the monomer densities Pfj, (r) as Gaussian func- 
tions, Pfj,{r) = S*^, exp[-7^(r - r^)^], taking for 5'^, 7^, 
and r^,^ — C, S, those values that provide the best fit 
of the simulation data, quoted in Fig. [2l Approximate 
expressions for the Pfi{k) are given in |22l |. The theo- 
retical results for $eff (fc) and hence $eff are shown in 
Fig. [4] along with the data extracted from the simulations. 
In view of the simplifying assumptions of Flory theory, 
the good qualitative agreement between simulations and 
theory is astonishing. The negative Fourier components 
[Figs.mjb) andUI^d)], are less pronounced in theory than 
in simulation, thus the former provides a lower threshold 
to the onset of clustering. 

The counterintuitive phenomenon of clustering in the 
complete absence of attraction might motivate experi- 
mental groups to assemble amphiphilic dendrimers in the 
lab. To this end let us summarize our guidelines for syn- 
thesizing clustering dendrimers. In a first step, suitable 
solvophobic core and solvophilic shell groups have to be 
chosen for the experiments, for which simulations on an 
isolated dendrimer are performed, leading to the core- 
and the shell-density profiles. While Flory theory pro- 
vides a reliable qualitative indicator whether the thresh- 
old to clustering has already been reached, full evidence 
can then be gathered by measuring the effective inter- 
actions in the more time-consuming simulations of two 
interacting dendrimers. From our observations for all 
dendrimers investigated, the following remarks should 
prove valuable: bigger end-groups and/or shorter end- 



group spacers lead to a stronger repulsion in $eff('') at 
short distances. Further, if the spacer length of the end- 
groups is increased and/or the distance between the two 
central particles and/or the end-group particle size is re- 
duced, <I'cff(?') becomes flatter or even develops a dip at 
small distances. The low dendrimer-generation number 
is encouraging for experimentalists because it hints at 
a rather straightforward synthesis process [l3|. Since 
^cE{r = 0) ^ k^T, clustering can easily be reahzed un- 
der ambient conditions in thermally activated processes. 

The findings of this work bear significance for soft mat- 
ter science and materials design at various levels. At the 
one-particle level, we have established that synthesizing 
open dendrimers with a segregated core-shell structure 
requires neither stiff bonds nor electrostatic repulsions 
as commonly believed: amphiphilicity is sufficient. At 
the many-body level, solutions of such dendrimers will 
display pronounced correlations at a single length scale, 
independently of the density (IQ] . allowing thus for well- 
controlled spatial modulation of confined liquids and, 
e.g., their local index of refraction, whose intensity can be 
tuned by changing the degree of confinement. Crystals 
formed by such systems show Aensiiy- independent lattice 
constants, a novel of self-assembly of condensed matter. 
Such crystals are quite unusual, since they are diffusive 
on the single-particle level, allowing thus mass transport 
but arrested, and thus rigid as a conventional solid, at the 
collective level [llj. Finally, on the fundamental level, we 
have demonstrated that within soft matter, bounded ef- 
fective interactions can be manipulated with the same 
degree of flexibility as diverging ones do. 

We thank D. Gottwald for helpful discussions. This 
work was funded by the FWF, Proj. No. 17823-N08 and 
by the DFG within the SFB-TR6, Project Section C3. 
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